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Continuous wavelet transform and cross-wavelet transform have been used to investigate the phase periodicity and syn- 
chrony of the monthly mean Wolf (Rz) and group (Rg) sunspot numbers during the period of June 1795 to December 
1995. The Schwabe cycle is the only one common period in Rg and Rz, but it is not well-defined in case of cycles 5-7 of 
Rg and in case of cycles 5 and 6 of _Rz. In fact, the Schwabe period is slightly different in Rg and Rz before cycle 12, but 
from cycle 12 onwards it is almost the same for the two time series. Asynchrony of the two time series is more obviously 
seen in cycles 5 and 6 than in the following cycles, and usually more obviously seen around the maximum time of a cycle 
than during the rest of the cycle. Rg is found to fit Rz better in both amplitudes and peak epoch during the minimum time 
time of a solar cycle than during the maximum time of the cycle, which should be caused by their different definition, and 
around the maximum time of a cycle, Rg is usually less than Rz . Asynchrony of Rg and Rz should somewhat agree with 
different sunspot cycle characteristics exhibited by themselves. 
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1 Introduction 

The single but most important index of solar activity has 
been the Zurich or Wolf sunspot number (the international 
sunspot number or the sunspot number) (Hathaway, Wil- 
son, & Reichmann 2002; Hathaway 2010). For more than 
100 years the Wolf sunspot number has served as the pri- 
mary time series to define solar activity (Hoyt & Schatten 
1998a, 1998b). It has been proven invaluable in studies of 
long-term variations in solar activity, especially as related to 
terrestrial climate (e.g., Eddy 1976; Hoyt & Schatten 1997; 
Li, Gao & Su 2005; ). 

The sunspot number is defined as i?^ = /s(10g + /), 
where k is the normalization constant for a particular ob- 
server, g is the number of sunspot groups, and / the number 
of individual sunspots visible over the solar disk. In spite 
of the apparent arbitrary nature of this formula, it has been 
found to correlate extremely well with other, more physi- 
cal measures of solar activity such as sunspot area, 10.7 cm 
radio flux, X-ray flare frequency, and magnetic flux (Hath- 
away, Wilson, & Reichmann 2002). Monthly values for 
are available from 1749 onward, but many of values 
between 1700 and 1850 are clearly based on inaccurate or 
missing data. Thus, these data are less accurate than more 
recent observations, and the R^ values appear to be too 
large by 25-50% prior to 1882 (Hoyt & Schatten 1995a, 
b, c, d; Wilson 1998; Faria et al. 2004). A new parame- 
ter, the group sunspot number (Rg) was introduced by Hoyt 
& Schatten (1998a, 1998b) as an alternative to the sunspot 



number. It uses only the number of sunspot groups but is 
normalized to make it to agree with the Zurich sunspot num- 
ber With this normalization the group sunspot number is 



given by Rg = ^^f^ where ki is the correc- 

tion factor for observer i, Qi is the number of sunspot groups 
observed by observer i, and N the number of observers 
used to form the daily value. The group sunspot number is 
more self-consistent and less noisy than the sunspot number. 
Daily, monthly, and annual observations were determined 
during the period of the years 1610- 1995. The Wolf and 
group sunspot numbers are the longest direct instrumental 
records of solar activity (Rigozo et al. 2001; Kane 2002; 
Ogurtsov et al. 2002; Echer et al. 2005). Hathaway, Wilson, 
& Reichmann (2002) examined the group sunspot number 
and compared the sunspot cycle characteristics it exhibits 
to those exhibited by the Wolf sunspot number. They found 
that the Wolf sunspot number continues to be valuable in 
capturing characteristics of the recent cycles that are not 
quite as well-reflected in the group number, and the group 
sunspot number is valuable for capturing the behavior in 
the earliest cycles that help to reveal long-term behavior 
Ogurtsov et al. (2002) used the yearly Wolf sunspot num- 
ber and the yearly group sunspot number to study long pe- 
riods of solar activity through Morlet wavelet analysis and 
Fourier analysis. Faria et al. (2004) compared the spectral 
features of the yearly group sunspot number and the yearly 
Wolf sunspot number with the use of multitaper analysis 
and Morlet wavelet analysis. Li et al. (2005) investigated the 
Schwabe and Gleissberg periods in the Wolf sunspot num- 
ber and the group sunspot number However, no one has in- 
vestigated phase relationship between them up to now, and 
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Fig. 1 The monthly mean group sunspot number (the black line) 
and the monthly mean sunspot number (the red line) from June 
1795 to December 1995. 



the difference between the two caused by their different def- 
inition has not yet been paid enough attention. 

According to the modern point of view, synchroniza- 
tion is an universal concept in nonHnear sciences (Pikovsky, 
Rosenblum & Kurths 2001; Maraun & Kurths 2004; Ro- 
mano et al. 2005). In the present study, we will use the cross- 
wavelet transform method to study phase relation between 
the monthly mean group sunspot number and the monthly 
mean sunspot number and investigate the influence of their 
different definition on the relation between them. 

2 Wavelet analyses of the monthly mean 
Wolf and group sunspot numbers 

2.1 Data 

The continuous time series of the monthly mean sunspot 
number (Rz) is available from January 1749, and the contin- 
uous time series of the monthly mean group sunspot number 
(Rg) is available from June 1795 to December 1995. Fig- 
ure 1 shows both Rg and R^ from June 1795 to December 
1995, which are both downloaded from the NOAA's web 
site'. These two kinds of data are used to investigate the 
phase relationship between them. As the figure shows, Rg 
obviously differs from R^ around the maximum time of a 
cycle, implying that they should be asynchronous. 

2.2 Continuous wavelet transform 
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Fig. 2 The continuous wavelet power spectra of the monthly 
mean group sunspot number (top) and the monthly mean sunspot 
number (bottom). The thick black contours indicate the 95% con- 
fidence level, and the white region below the thin solid line indi- 
cates the cone of influence (COI) where edge effects might distort 
the picture (for details, see Torrence & Compo 1998). 



age for detecting the localized and (pseudo-) periodic fluctu- 
ations by using the limited time span of the data (Torrence 
& Compo 1998; Li et al. 2005). Here the Morlet continu- 
ous wavelet transform (CWT) is used with its dimensionless 
frequency ojq = 6. Figure 2 shows the continuous wavelet 
power spectra of Rg and R^. There are evidently common 
features in the wavelet power spectra of the two time se- 
ries. The both have a large-scale periodicity (about 1 1 years, 
namely the Schwabe cycle) of the highest power, above the 
95% confidence level. However, the Schwabe cycle does not 
significantly exist in cycles 5 to 7 for Rg and in cycles 5 and 
6fori?2. 

Figure 3 displays the period length of the Schwabe cycle 
varying with time respectively for the two time series. At a 
certain time point, the Schwabe period (scale) of Rg (or R^) 
has the highest spectral power among all considered time 
scales in its local wavelet power spectrum, thus the Schwabe 
period length is determined at a certain time point. As the 
figure shows, the Schwabe period for Rg and R^ actually 
differs from each other before cycle 12, but from cycle 12 
onwards it is almost the same for the two time series. If 
period (frequency) of two time series is different, the two are 
asynchronous. The differences of the Schwabe cycle length 
for the two time series should lead to a phase asynchrony 
between them. 



Wavelet analysis involves a transform from an one-dimensionali 

time-series to a diffuse two-dimensional time-frequency im- ""^ Cross-wavelet transform 



' ftp - // ftp.ngdc.noaa.gov/STP/SOLAR-DATA/ 
SUNSPOT- NUMBERS 



Cross-wavelet transform (XWT) is an extension of wavelet 
transforms to expose the common power and relative phase 
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Fig. 3 The period length of the Schwabe cycle for the monthly 
mean group sunspot number (the red line) and that for the monthly 
mean sunspot number (the black solid line). The two dashed lines 
indicate the COL 



between two time series in time-frequency space (Li et al 
2009 and references therein). The cross - wavelet transform 
of two time series X and Y is defined as W^^ — W^W^*M 
where * denotes complex conjugation and W'^ and are 
the continuous wavelet transforms of the individual time se- 
ries (Grinsted et al. 2004). The complex argument arg{W'^^)\ 
can be interpreted as a local relative phase between X and Y 
in time-frequency space, namely the phase angle difference 
of X and Y. We employ the codes provided by Grinsted 
et al. (2004)^ to get the XWT of the two time series in or- 
der to investigate their phase relationship. Figure 4 shows 
the XWT of the two time series. The two time series are in 
phase in the area, namely about an 11 yr periodic belt, with 
significant common power, because all arrows point right in 
the area. 

The real and imaginary parts of a wavelet power spec- 
trum are considered separately, and the phase angle at any 
point of the wavelet power spectrum can be obtained through^ 
its complex spectral amplitude (Li et al. 2009). Figure 5 
shows the phase angle (argument) of the complex cross- 
wavelet amplitudes varying with frequency and time. It is 
the relative phase of the two time series. The results shown 
in the figure clearly indicate that the relative phases coher- 
ently fluctuate within a rather large range of frequencies, 
which correspond to the period scales of larger than 4 years. 
At smaller scales, there are no regular oscillatory patterns in 
the corresponding frequency bands, and the high-frequency 
components demonstrate a noisy behavior with strong phase 
mixing. 



^ http://www.pol.ac.uk/home/research/waveletcoherence/ 
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Fig. 4 The cross wavelet transform of the monthly mean group 
sunspot number and the monthly mean sunspot number. The thick 
black contours indicate the 95% confidence level, and the region 
below the thin line indicates the COL The relative phase relation- 
ship is shown as arrows with in-phase pointing right, anti-phase 
pointing left, and the former leading the latter by 90° pointing 
straight down. 
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Fig. 5 The relative phase angles (arguments) of the complex 
cross-wavelet amplitudes varying with frequency and time. The 
thin solid line indicates the COL 
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Fig. 6 Tiie average (tire thick solid line) of relative phases over 
all period scales shown in Figure 5 varying with time. The dashed 
(solid) vertical lines show the minimum (maximum) times of cy- 
cles. 



Based on Figure 5, we calculate the averages of relative 
phase angles over all scales for all time points of the consid- 
ered interval, which are displayed in Figure 6. Their corre- 
sponding standard deviations are also calculated and shown 
in Figure 7. As the two figures show, the relative phase of 
the two time series fluctuates within larger amplitudes in cy- 
cles 5 and 6 than that in all the following cycles does, and it 
fluctuates within larger amplitudes usually around the max- 
imum time of a cycle than around the minimum time of the 
cycle. Asynchrony of the two time series is thus more ob- 
viously seen in cycles 5 and 6 than in the following cycles, 
and usually more obviously seen around the maximum time 
of a cycle than during the rest of the cycle. 

For the two time series, the mean relative phases are 
shown in Figure 8 in dependence on the considered fre- 
quencies. The corresponding standard deviations have been 
calculated as well and are given in the second panel of the 
figure. As the figure shows, the relative phases obviously co- 
herently fluctuate within time scales of about 8 to 13 years 
(and even within 8 to 25 years), where their relative phases 
are within ±0.2 radian, and their standard deviations are 
less than 0.5 radian. Although the relative phases are always 
negative within time scales of about 8 to 13 years, seemingly 
indicating leading Rg, the corresponding standard devi- 
ations are so large that we could not think R^ leading Rg . At 
time scales of less than 8 years, the relative phase angles vi- 
olently fluctuate and their standard deviations are very large. 
The relative phase angles increase after time scales of larger 
than 25 years, and their standard deviations also increase. 

Taking very small or very large reference time scales 
into account should lead to non-coherent behavior of the 
phase variables assigned to Rg and R^. In the following. 



ffl 1.5 



.9 1 
■> 

0) 

°0.5 




1800 1810 1820 1830 1840 1850 1860 1870 1880 1890 
Calendar Year 




1900 



1920 1940 1960 
Calendar Year 



1980 



2000 



Fig. 7 The corresponding standard deviation (the thick line) of 
the average of relative phases shown in Figure 6. The dashed 
(solid) vertical lines show the minimum (maximum) times of cy- 
cles. 
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Fig. 8 Angle difference (the relative phase) of the XWT power 
spectrum of of the monthly mean group sunspot number and 
the monthly mean sunspot number as a function of periods (the 
top panel) and its corresponding standard deviations (the bottom 
panel). Positive values should be interpreted as the monthly mean 
group sunspot number leading the monthly mean sunspot number. 
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Fig. 9 The average of relative phases over the period scales re- 
spectively of 8 to 13 years (the solid line) and of 8 to 25 years (the 
dotted line) varying with time (the top panel), and its correspond- 
ing standard deviation (the lower panel) with the XWT method 
used. Solid (dotted) line corresponds to solid (dotted) line in the 
two panels. 
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Fig. 10 Relationship between Rg and R- . Rg is divided into two 
group: in one (Gl) -Rg is less than its median (the yellow asterisks), 
and in the other (G2) Rg is larger than the median (the blue aster- 
isks). A linear fit is done for each of the two groups (the dashed 
line is for Gl, and the solid for G2). The dotted line is the line of 
that Rg = R^. 



we therefore focus attentions on periodicities in the coher- 
ent range around the Schwabe cycle as the reference time 
scales. Because the high-frequency components of Rg and 
Rz are filtered out in such a case, the resulting phases are 
well-defined and may be used to study the varying relation- 
ship between Rg and Rz- Based on Figure 5, we calculate 
the averages of relative phase angles over scales of 8 to 13 
years for all time points of the considered interval, which are 
displayed in Figure 9. Their corresponding standard devia- 
tions are also calculated and shown in the figure. As the fig- 
ure indicates, for all cycles except cycles 5 and 6, the mean 
relative phase angles are very small (within ±0.35 radian), 
fluctuating around the value of zero, and the standard devia- 
tions are small (less than 0.5 radian). It is the low-frequency 
components of hemispheric flare activity in period scales 
around the Schwabe cycle that are responsible for the strong 
phase synchronization. Based on Figure 5, we calculate the 
averages of relative phase angles over scales of 8 to 25 years 
for all time points of the considered interval, which are also 
displayed in Figure 9. Their corresponding standard devia- 
tions are calculated and shown in the figure. As the figure 
indicates, for all cycles except cycles 5 and 6, the mean rel- 
ative phase angles are very small (even within ±0.2 radian), 
fluctuating around the value of zero, and the standard devi- 
ations are small (less than 0.7 radian). For cycles 5 and 6, 
there are no coherent frequencies. 



3 Influence of their different definition on the 
relation between them 

As we know, during the minimum time of a cycle a group 
of sunspots has only one sunspot in the most cases or two 
sometimes, thus (lOg + f)/g « 11. However, during the 
maximum time of a cycle a group of sunspots has one to 
two sunspots sometimes, or several (or even more) sunspots 
sometimes, thus (lOg + is larger than 1 1 and believed 
even usually larger than 12.08 (the normalization constant 
in the definition of Rg). Therefore, there should certainly 
exist difference between the two caused by their different 
definition. The median value of Rg is about 46.4, and Rg 
is divided into two group: in one (Gl) Rg is less than, and 
in the other (G2) Rg is larger than the median. A linear fit 
is done for each of the two groups, which is shown in Fig- 
ure 10. The fitting line is 7?g = 1.5 + 0.897 x R^ with 
the coiTelation coefficient of 0.965 for the group Gl, and 
Rg = 8.2 + 0.883 X Rz with the correlation coefficient of 
0.962 for the group G2. As the figure shows, the fitting line 
is very close to the straight line Rg = Rz for Gl. How- 
ever, for G2 the fitting line is generally below the straight 
line Rg = Rz, and the deviation between the two increases 
with the increase of Rg . Conclusively, Rg matches Rz very 
well around the minimum time of a cycle, but during the 
maximum time of a cycle it seems usually less than Rz- Rg 
and Rz are doomed by their different definition to present 
different relationship between them at different stage of a 
cycle. 
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Fig. 11 Tlie top panel: difference (asterisks) between the max- 
imum time of the sunspot number and the group sunspot num- 
ber in a cycle. The dashed line shows no differences. The bottom 
panel: the maximum values of the sunspot number (asterisks) and 
the group sunspot number (crosses). The solid line shows the av- 
erage (135.5) of the sunspot number maxima over cycles 10 to 22, 
and the dashed, the average (132.0) of the group sunspot number 
maxima. 



The maximum value of a sunspot cycle is an important 
parameter to characterize the sunspot cycle. Figure 1 1 dis- 
plays the maximum values of the Wolf sunspot number and 
the group sunspot number in the modern cycles (cycles 10 to 
22). The maximum value of the sunspot number in a cycle is 
usually larger than that of the group sunspot number (cycles 
10 to 12 and 18 to 22) or very close to the corresponding one 
of the group sunspot number (cycles 13 to 17), confirming 
the aforementioned result. The average of the sunspot num- 
ber maxima over cycles 10 to 22 is 135.5, and the average 
of the group sunspot number maxima is 132.0, slightly less 
than the former. We also show in the figure the difference 
between the maximum time of the sunspot number and the 
group sunspot number in a cycle for the modern cycles. In 
9 of the total 13 cycles, the difference is obvious, indicating 
asynchrony of Rg and . 

4 Conclusions and discussions 

In the present study, the continuous wavelet transform and 
the cross-wavelet transform have been proposed to inves- 
tigate the phase synchrony of the monthly mean Wolf and 



it is 10.7 years. However, in cycles 5 to 7 the Schwabe cycle 
does not significantly exist for Rg, and in cycles 5 and 6 for 
Rz- 

The Schwabe period for Rg and R^ actually slightly dif- 
fers from each other before cycle 12, but from cycle 12 on- 
wards it is almost the same for the two time series. If pe- 
riod (frequency) of two time series is different, the two are 
asynchronous. The slight differences of the Schwabe cycle 
length for the two time series should lead to a slight phase 
asynchrony between them. 

The cross-wavelet transform (XWT) of the two time se- 
ries shows that there is an area, locating at the Schwabe- 
cycle periodic belt, where the two time series are approx- 
imately in phase. Through the XWT analysis, the relative 
phase angles are found to coherently fluctuate in a small 
angle range only within a range of frequencies which corre- 
sponds to time scales around the Schwabe cycle, and even 
within time scales of about 8 to 25 years. The high-frequency| 
components demonstrate a noisy behavior with strong phase 
mixing. It is the low-frequency components of Rg and Rz 
in period scales around the Schwabe cycle that are responsi- 
ble for the strong phase synchronization. Taking very small 
or very large reference time scales into account should lead 
to non-coherent behavior of the phase variables assigned to 
Rg and R^ at the respective scales. Asynchrony of the two 
time series is found more obvious in cycles 5 and 6 than in 
the following cycles, and more obvious usually around the 
maximum time of a cycle than during the minimum time 
of the cycle, indicating that Rg should fit R^ better in both 
amplitudes and peak epoch during the minimum time time 
of a solar cycle than during the maximum time of the cycle. 
As we know, during the minimum time of a cycle a group of 
sunspots has only one sunspot generally or two sometimes, 
thus {Wg + f) / g ~ 11. However, during the maximum 
time of a cycle a group of sunspots has one to two sunspots 
sometimes, or several (or even more) sunspots sometimes, 
thus (lOg + /) /g is larger than 1 1 and believed even usually 
larger than 12.08 (the normalization constant in the defini- 
tion of Rg), and further it is scattered within a wider range 
than during the minimum of the cycle. Rg is destined to 
match Rz better during the minimum time than during the 
maximum time of a cycle, making the two more obviously 
asynchronous during the maximum time than during the 
minimum time of a cycle, and further, Rg is usually less 
than Rz around the maximum time of a cycle. 

Somewhat in agreement with the more obvious asyn- 
chrony of the two during the maximum times than during 
the the corresponding minimum times of solar cycles, the 
"Waldmeier Effect" - the anti-correlation between cycle am-f 
plitude and the elapsed time between minimum and maxi- 



group sunspot numbers, the longest direct instrumental records|wMm of a cycle - is much more apparent in the Zurich num- 



of solar activity, and the data used here span from June 1795 
to December 1995. 

The continuous wavelet transforms of the two time se- 
ries show that the 11 -year Schwabe cycle is the only one 
period of statistical significance for the two time series, and 



bers; the "Amplitude-Period Effect" - the anti-correlation 
between cycle amplitude and the length of the previous cy- 
cle from minimum to minimum - is also much more apparent 
in the Zurich numbers; the "Even-Odd Effect" - in which 
odd-numbered cycles are larger than their even-numbered 
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precursors - is somewhat stronger in the Group numbers References 
but with a tighter relationship in the Zurich numbers; the 
'Secular Trend' - the increase in cycle amplitudes since the 
Maunder Minimum - is much stronger in Group numbers 
(Hathaway, Wilson & Reichmann 2002). Due to that Rg 
matches worse during the maximum time than during 
the minimum time of a cycle, they exhibit the above dif- 
ferent sunspot cycle characteristics and more obvious asyn- 
chrony during the maximum time of the cycle. 

Sunspot number, group sunspot number and sunspot area§ 
are widely and frequently utilized in astronomy and Earth 
science to embody long-term variations of solar activity. 

However, reviewed from the phase relationship between each| Hoyt D V. Schatten K H 
two of the three, they are slightly asynchronous in phase, al- 
though they are coherent in the low-frequency components 
corresponding to the period scales around the Schwabe cy- 
cle, this is inferred the main reason why "the Zurich num- 
bers follow the 10.7-cm radio flux and total sunspot area 
measurements only slightly better than the Group numbers" 
(Hathaway, Wilson & Reichmann 2002). Which one of the 
three can be served as the best primary time series to define 
solar activity? it is an open question, and careful attention 
should be paid. 
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